REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  NO.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments 

regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 

Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington  VA,  22202-4302. 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  oenalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 
14-12-2009 


4.  TITLE  AND  SUBTITLE 

Controlled  dynamic  fragmentation  of  ceramics 

Final  Technical  Report 


2.  REPORT  TYPE 
Final  Report 


3.  DATES  COVERED  (From  -  To) 

18-M-2005  -  17 -Nov-2 007 


5a.  CONTRACT  NUMBER 
W91  INF-05-1-0370 


5b.  GRANT  NUMBER 


6.  AUTHORS 
Prof.  J.F.  Molinari 


5c.  PROGRAM  ELEMENT  NUMBER 
611102 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRESSES 

Johns  Hopkins  University 
Johns  Hopkins  University 
3400  N.  Charles  St. 

Baltimore,  MD  21218  -2686 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND 
ADDRESS(ES) 

U.S.  Army  Research  Office 
P.O.Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 
ARO 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

48389-MS. 1 


12.  DISTRIBUTION  AVAILIBILITY  STATEMENT 
Approved  for  Public  Release;  Distribution  Unlimited 


13.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  contrued  as  an  official  Department 
of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


14.  ABSTRACT 

Fragmentation  is  a  key  damage  mechanism  determining  ballistic  impact  performance  of  ceramic  armors  as  well  as 
reliability  of  gun  barrels.  This  program,  which  was  sponsored  by  the  Army  Research  Office,  focuses  on  the  relation 
between  defect  populations  and  fragment-size  distributions  via  robust  physics-based  numerical  simulations.  The 
originality  of  the  present  simulations  lies  in  the  capability  of  dealing  with  fragmentation  of  large  and  heterogeneous 
structures.  The  understanding  of  how  the  microstructure  acts  on  the  global  macroscopic  response  is  our  main 


15.  SUBJECT  TERMS 

Dynamic  fragmentation;  defects;  cohesive  elements; 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

15.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF  PAGES 

Jean-Francois  Molinari 

UU 

UU 

UU 

UU 

19b.  TELEPHONE  NUMBER 

412-169-3241 

Standard  Form  298  (Rev  8/98) 
Prescribed  by  ANSI  Std.  Z39. 1 8 


REPORT  DOCUMENTATION  PAGE  (SF-298  W911  NF-08-1-0150) 
FINAL  TECHNICAL  REPORT 


Controlled  dynamic  fragmentation  of  ceramics 
PI:  JF  Molinari  (Johns  Hopkins  University,  Ecole  Polytechnique  Federate  of  Lausanne) 


Attention:  Dr.  D.  Stepp,  Army  Research  Office 


Fragmentation  is  a  key  damage  mechanism  determining  ballistic  impact  performance  of  ceramic  armors  as  well 
as  reliability  of  gun  barrels.  This  program,  which  was  sponsored  by  the  Army  Research  Office,  focuses  on  the 
relation  between  defect  populations  and  fragment-size  distributions  via  robust  physics-based  numerical 
simulations.  The  originality  of  the  present  simulations  lies  in  the  capability  of  dealing  with  fragmentation  of 
large  and  heterogeneous  structures.  The  understanding  of  how  the  microstructure  acts  on  the  global  macroscopic 
response  is  our  main  concern.  Defects  are  strategically  placed  during  the  design  of  a  material  or  a  structure  in 
order  to  optimize  the  geometrical  (shape  and  size)  or  velocity  characteristics  of  fragments. 

Key  simulation  features: 

•  Comprehensive  physics  to  capture  dynamic  crack  propagation 

•  Presence  of  material  defects 

•  Novel  modeling  strategies:  Discontinuous  Galerkin  fragmentation 

•  Massively  parallel  simulations 


Below,  we  detail  a  few  salient  scientific  results  of  the  program.  The  complete  list  of  publications  is  given  at  the 
end  of  the  document. 


Part  1:  Fragment  sizes  resulting  from  brittle  fragmentation. 

The  average  fragment  mass  is  perhaps  the  most  studied  physical  property  resulting  from  dynamic  fragmentation. 
This  is  no  surprise,  as  fragment  masses  find  applications  in  several  engineering  and  scientific  fields  (including 
ballistic  impact  performance).  Starting  from  the  pioneering  analytical  work  of  Mott  [1],  up  to  the  seminal  energy 
model  of  Grady  [2],  several  models  have  sought  to  understand  the  link  between  material  properties,  strain  rate, 
and  fragment  mass.  Recent  work  has  focused  on  resorting  to  numerical  simulations  to  address  the  limitations  of 
analytical  work.  Our  contribution  follows  this  trend.  Our  originality  lies  in  coupling  the  best  practices  in 
mechanical  modelling  of  rupture  with  a  rigorous  analysis  of  the  physics  at  play.  In  our  approach,  material  rupture 
is  captured  by  having  recourse  to  cohesive  elements  [3],  which  are  dynamically  inserted  in  a  previously 
continuous  finite-element  mesh.  The  influence  of  material  defects  having  been  seldom  studied,  we  conducted 
parametric  studies  of  fracture  parameters  distributions  to  capture  statistical  heterogeneity. 

Our  results  reveal  that,  despite  statistical  heterogeneity  at  the  small  scale,  and  complex  dynamics,  the  output  of 
fragmentation  is  not  random  but  displays  universal  features,  see  figure  1.  The  normalisation  procedure  that  we 
used  to  obtain  the  remarkably  simple  law  is  the  subject  of  a  publication  [4].  It  involves  three  parameters,  which 
characterize  the  distribution  of  defects.  They  are  the  largest  defects  (also  known  as  the  weakest  links),  the 
average  defect  size,  and  the  left  slope  of  the  defect  distribution.  We  have  demonstrated  that  this  slope  controls 
the  rate  of  crack  initiation  at  a  given  strain  rate.  These  parameters  being  known,  we  are  able  to  predict  the 
average  fragment  mass  obtained  from  brittle  fragmentation  events.  The  inverse  problem  is  perhaps  even  more 
interesting:  our  theoretical  model  opens  the  door,  from  inspecting  the  average  fragment  mass  at  various  strain 
rates,  to  characterizing  the  defect  distribution  in  a  given  material. 

In  a  follow  up  study  [5],  we  have  analyzed  the  fragment  size  distributions.  The  end  result  of  a  brittle 
fragmentation  event  is  a  diverse  population  of  fragment  sizes.  We  are  for  instance  well  acquainted  with  the 
breakage  of  kitchen  plates,  resulting  in  large  pieces,  which  are  easily  cleaned,  and  finer  wreckage,  which  may  go 
unnoticed  for  weeks.  Several  semi-empirical  models  describing  this  diversity  in  fragment  masses  have  been 
derived  in  the  last  century.  We  compared  our  numerical  predictions  to  these  models,  and  discovered  a  surprising 
feature.  If  we  normalize  the  fragment  sizes  with  the  average  fragment  size  (which  we  can  now  reliably  predict, 
see  above),  all  our  distributions  fall  under  a  unique  mixed-power-exponential  form  (fig.  2a).  The  end  result  of 


REPORT  DOCUMENTATION  PAGE  (SF-298  W911  NF-08-1-0150) 
FINAL  TECHNICAL  REPORT 


this  numerical  investigation  is  that  we  are  not  only  in  position  to  predict  an  average  fragment  mass,  but  also  the 
statistics  of  the  entire  defect  distribution. 

Interestingly,  the  largest  fragments  do  not  seem  to  obey  such  a  simple  scaling,  as  we  see  in  fig.  2b.  The  influence 
of  defects  on  the  emergence  of  large  fragments  remains  a  mystery  that  we  are  currently  studying  through 
extreme  value  statistics  analysis. 
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Fig  1:  Schematic  of  ceramic  ring  expansion  test  showing  the  initiation  of  multiple  fracture  sites  and  the 
development  of  a  complex  network  of  stress  waves  interactions.  On  the  right,  our  results  show  the  surprising 
robustness  of  the  average  fragment  size  prediction  as  a  function  of  strain  rate  for  a  ceramic  material.  The 
normalization  procedure  of  the  axes  is  subject  of  a  publication,  and  involves  an  explicit  dependence  on  materials 
properties  (elastic  and  fracture)  and  defect  population  statistics.  It  is  noteworthy  that  all  the  defect  distributions 
that  were  considered  (Weibull,  Gauss  and  Uniform)  fall  on  a  single  universal  curve. 


Fig  2  :  a)  Inverse  cumulative  density  function  of  fragment  size  (s).  Each  test  involves  different  material 
properties,  defect  distribution,  and  strain  rates.  Yet,  the  normalisation  by  the  average  fragment  size  yields  a 
unique  master  curve,  b)  Extreme  value  statistics  reveal  that  the  same  simple  scaling  law  does  not  hold  true 
for  the  largest  fragments.  The  link  between  defect  population  and  the  largest  fragments  remains  a  mystery’. 


Part  2:  Discontinuous  Galerkin  fragmentation  and  massively  parallel  simulations 


In  collaboration  with  Prof.  R.  Radovitzky  and  A.  Seagraves  (MIT)  we  have  conducted  parallel  simulations  of  3D 
fragmentation  (see  figure  3)  using  a  Discontinuous  Galerkin  (DG)  approach  [7],  which  was  developed  at  MIT 
[6,8]. 

The  DG  approach  is  finding  an  increasingly  broader  appeal  in  the  mechanics  community,  due  to  a  more  accurate 
description  of  shock  fronts,  than  what  is  achieved  with  regular  continuous  Galerkin  approaches  (regular  finite 
element  codes).  In  DG,  all  element  boundaries  are  initially  separated,  and  flux  terms  are  added  to  satisfy 
equilibrium  weakly.  This  initial  splitting  of  all  mesh  elements  fixes  the  mesh  topology.  In  standard  cohesive 
elements  codes,  insertion  of  cracks  must  be  done  dynamically,  which  implies  that  the  mesh  topology  is 
constantly  changing.  This  is  easily  achieved  on  single  processor  machines,  but  is  a  much  harder  endeavour  on 
multiple  processor  machines,  in  which  cracks  may  propagate  from  a  processor  to  another. 
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DG  methods  open  the  door  to  massively  (and  scalable)  fragmentation  simulations  of  3D  structures. 

Therefore,  the  method  shows  excellent  potential  for  scaling  on  HPC  machines.  The  PI  has  purchased  significant 
computing  power  and  DG  is  currently  tested  to  generated  fragment  size  statistics  in  simulations  involving 
hundreds  of  processors  (see  figure  3).  We  are  currently  evaluating  how  defects  influence  fragment  shapes. 
Simulations  reveal,  as  fragmentation  damage  unfolds,  a  complex  yet  “organized”  percolation  network  of  cracks. 
Strategically  placed  defects  influence  the  organization  of  this  network  and  can  effectively  control  the  sharpness 
of  fragments  in  impacted  ceramics. 


Figure  3:  a)  3D  ceramic  plate  under  rapid  biaxial  tension  (104/s);  a  random  Weibull  distribution  of  initial 
defects  yields  heterogeneous  fragment  sizes  or  roughly  equi-axed  shapes,  b)  Parallel  simulations  of  a 
fragmenting  concrete  dome  due  to  an  internal  explosion.  Each  colour  represents  a  different  processor.  Notice 
how  cracks  propagate  from  one  processor  to  another. 


Key  achievements: 

Robust,  yet  simple,  numerical  and  analytical 
predictors  of  fragments  sizes,  which  are  function  of 
material  parameters  and  defect  distributions. 

Ability  to  conduct  massively  parallel  fragmentation 
simulations. 
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